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Abstract 

The non linear sigma model O(N) (g) 0(2) /0(N — 2)<S> 0(2) describing the phase 
transition of N-components helimagnets is built and studied up to two loop order 
in D = 2 + e dimensions. It is shown that a stable fixed point exists as soon 
as N is greater than 3 (or equal) in the neighborhood of two dimensions. The 
critical exponents v and r] are obtained. In the N = 3 case, the symmetry of 
the system is dynamically enlarged at the fixed point from 0(3) (g> 0(2)/0(2) to 
0(3) (g> 0(3)/ 0(3) ~ 0(4)/0(3). We show that the order parameter for Heisenberg 
helimagnets involves a tensor representation of 0(4) and we verify it explicitly at 
one loop order on the value of the exponents. We show that for large N and in 
the neighborhood of two dimensions this nonlinear sigma model describes the same 
critical theory as the Landau-Ginzburg linear theory. As a consequence, the critical 
behaviour evolves smoothly between D = 2 and D = 4 dimensions in this limit. 
However taking into account the old results from the D = 4 — e expansion of the 
linear theory, we show that most likely the nature of the transition must change 
between D = 2 and D = 4 dimensions for small enough ./V (including N = 3). The 
simplest possibility is that there exists a dividing line N C (D) in the plane (N,D) 
separating a first-order region containing the Heisenberg point at D = 4 and a 
second-order region containing the whole D = 2 axis. We conclude that the phase 
transition of Heisenberg helimagnets in dimension 3 is either first order or second 
order with 0(4) exponents involving a tensor representation or tricritical with mean 
field exponents. 
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1 Introduction 



In the recent past, the critical behavior of frustrated spin systems has been the subject of 
intensive theoretical, numerical and experimental studies (see [?][?] and references therein 
and [?]). However, there is still no definite conclusion about the nature of the phase tran- 
sition that occurs in these systems. One of the most striking features of frustrated models 
is their non trivial ground state which for continuous spin models is in general a canted 
ground state. Well known examples are incommensurate helimagnets and triangular anti- 
ferromagnets among others. As a consequence of the non collinear ordering the 0(3) spin 
rotation group is completely broken in the low temperature phase so that the relevant 
order parameter is a rotation matrix instead of a vector as in ferromagnetic-like models. 
One may thus wonder if canted spin models belong to a new universality class. Up to 
now, no definite answer is known. Experiments on rare earth helimagnets such as Ho, 
Tb, or Dy for example, do not show any clear evidence for a universal critical behaviour. 

From a theoretical point of view, early renormalization group (RG) studies by Garel 
and Pfeutyf?] found no stable fixed point in the neighborhood of D = 4 by studying 
the Landau-Ginzburg theory of a commensurate helimagnet. This is an example of a 
fluctuation-induced first order transition. However, early Monte-Carlo studies of several 
canted models [?][?] found evidence for a continuous transition in three dimensions. If 
taken seriously these numerical results are in contradiction with the 4 — e prediction. Of 
course, it is notoriously difficult to discriminate between first and second order phase 
transitions in Monte-Carlo studies (the two-dimensional five-states Potts model is a well- 
known weird case for example). Strictly speaking, one cannot exclude the existence of a 
stable fixed point which manifests itself at a finite distance from D — 4, unreachable in 
an e expansion around D — 4. If this happens to be true, then standard perturbative 
methods are of no use to study the critical behavior of canted spin models. 

There is an alternative perturbative approach to this problem which is the low tem- 
perature expansion of the non-linear sigma (NLa) model suited to the symmetry breaking 
scheme of these canted models. In this paper we focus on a simple commensurate heli- 
magnet which is the triangular antiferromagnet with N-component classical spins. By 
stacking triangular planes, this magnet exists in all integer dimensions (D > 2). 

In the case of the triangular antiferromagnet (AFT) with Heisenberg classical spins 
i.e. N — 3, the massless modes live in the homogeneous space G/H = 0(3) ® 0(2) /O (2). 
Some results from the D = 2 + e expansion of a NLa model based on this coset G/H 
have been recently reported[?] [?]. It has been found that up to two loop a stable fixed 
point which is the N = 4 Wilson-Fisher fixed point shows up in the vicinity of D — 2. 
Thus no new universality class is required in the case of canted spin models. One meets 
the general phenomenon of increased symmetry at a critical point since at this point the 
model is 0(3) <g> 0(3) = 0(4) instead of 0(3) <g> 0(2) symmetric. 

In this paper, we extend our previous analysis to the case with N > 3 components. We 
build up the relevant nonlinear sigma model and analyze its RG properties by standard 
field-theoretic techniques. 

If one believes that both the e = 4 — D and e = D — 2 perturbative results can 
be extented to non-zero e, in the neighborhood of D = 2 and D = 4, as it is the case 
for the 0(N) models, the simplest hypothesis which agrees with both e expansions is 
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the following: there is a tricritical surface separating the basin of attraction of the 0(4) 
fixed point found near D = 2 from a first order runaway region found in the vicinity 
of D = 4. The phase transition of canted magnets is thus either first order or second 
order with 0(4) or tricritical exponents (i.e. mean-field in D = 3). This hypothesis has 
been previously proposed in ref.[?]. Recent extensive Monte-Carlo studies [?] performed 
directly in D=3 point towards a second-order transition with 0(4) exponents, a fact that 
was missed by previous lower-statistics studies. This means that presumably the critical 
surface for N = 3 lies between D=3 and D=4. 

The manifold G/H is topologically equivalent to 0(3) but as metric spaces they are 
different. The RG properties of the corresponding non-linear sigma model are a priori 
sensitive to the metric properties. However the study of purely topological properties can 
be performed directly on 0(3) as in ref.[?]. The study of defects reveals the presence 
of Z 2 vortices that are probably liberated in the high-temperature phase of the strictly 
two-dimensional AFT model. In this work we will ignore global aspects and concentrate 
on configurations with zero vorticity, leaving for the future the study of the defects on 
the phase transition. 

In this paper, we present the detailed renormalization group study of canted spin 
systems in D = 2 + e. In section I we show how the effective continuum action is obtained 
from a lattice Hamiltonian with Heisenberg spins. In section II the group theoretical 
construction of the non-linear sigma (NLcr) model is presented. In section III the two 
loop recursion relations as well as the Callan-Symanzik 7-function are given. The critical 
exponents v and rj are calculated. Special attention is given to the nature of the order 
parameter which is shown to belong to the tensor representation of 0(4). In section IV 
known results from both 4 — e and 1/N expansions are recalled for convenience. They are 
discussed and compared with the 2 + e results in section V. Our conclusions are contained 
in section VI. 

2 Continuum limit and effective action 
2.1 General analysis 

The effective action that describes the long distance behavior of a lattice model is obtained 
by taking the continuum limit of the microscopic Hamiltonian: 

H = -^JijSi.Sj. (1) 

In this equation the vectors Sj are classical Heisenberg spins with fixed unit length. In 
a ferromagnetic system, this continuum limit is achieved by letting the spins S fluctuate 
around their common expectation value. Relative fluctuations between neighboring spins 
are assumed to be smooth enough so that we may replace Sj.Sj by (VS(x)) 2 . 

When the interaction distribution { J^} leads to a canted ground state the continuum 
limit is less obvious since neighboring spins do not fluctuate around the same mean ex- 
pectation value. To overcome this difficulty, one has to consider the magnetic cell with 
n sublattices (S 1 , . . . , S n ) as the basis of a new superlattice where the continuum limit 
is taken. Practically, this procedure depends on the detailed microscopic model: lattice 
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symmetry, ground state structure and interaction parameters. We shall however give 

qualitative arguments valid for many canted models. 

Let us define in each elementary cell an orthonormal basis {e a (x)}: 

e a (x).e b (x) = 5 ab ; a= 1,2,3 , (2) 

where x is a superlattice index. We may parametrize our n sublattice spins S a (x),a = 
1, .., n as: 

S a (x)=Y,C:(x)e a (x). (3) 



In the ground state, all the S Q are in general not independent. There is in fact a maximum 
of three of them which are independent. Equivalently, there is a minimum of n — 3 
linear combinations of the S a (x) which have zero expectation value in the ground state. 
Such combinations cannot be part of an order parameter. They correspond to relative 
motions of the spins within each unit cell. They are massive modes with short range 
correlations and are thus irrelevant to the critical behavior. We ignore them by imposing 
the constraints that locally, i.e. within each unit cell, the spins are in the ground state 
configuration. We call this requirement "local rigidity" . Thus, up to an appropriate field 
redefinition, the order parameter of canted magnets will be the orthonormal basis {e a (x)} 
defined on each site of the superlattice. As a consequence canted magnets are equivalent 
in the critical region to a system of interacting solid rigid bodies. The continuum effective 
action Si may be obtained through the standard gradient expansion of the e a (x) as in 
ferromagnets: 

Si = l jd D x f> a (Ve a ) 2 , 

1 J a=i (4) 



e a (x).e b (x) = 5 ab , 
where the ground state is given by the minimization equations: 

Ve a (x) = 



(5) 



e a (x) = e 

The p a , a = 1, 2, 3 are coupling constants which depend on the particular lattice model 
we started with. The partition function Z is : 

Z = J De lj2i3 (x) (j[ 5(e a (x).e b (x) ~ <U)) e~ s ^ T . (6) 

\ ab / 

When two couplings p a vanish, one recovers, integrating over the corresponding e a , the 
action of the standard non linear sigma model 0(3)/0(2) corresponding to collinear ferro- 
or anti- ferromagnets. In all the other cases, among the nine fields e a (x), taking into 
account the constraints (|2|), one sees that there are three independent fluctuating fields 
corresponding to the three Goldstone modes, or spin waves, resulting from the breakdown 
of the 0(3) group. Each one corresponds to infinitesimal rotations around each of the 
e a (x)'s. The couplings p a are the associated stiffness constants which depend on the 
detailed microscopic model. They are deeply connected to the symmetry properties of 
the lattice Hamiltonian as we shall see. 
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In addition to the usual 0(3) rotational invariance, the symmetry group G of the 
Hamiltonian contains, in general, a discrete group T of transformations {T s } mixing 
together the sublattices S Q , or equivalently the e a . These transformations may belong to 
the space group of the lattice, as in triangular antiferromagnets, but may be some more 
complicated objects, such as "gauge" transformations as in the Villain latticef?]. The 
order parameter {e a }, a = 1, 2, 3 thus transforms under G as: 

<:,-Y,c,,<!,: £/gO(3), 

3 

(7) 

e a - ]T( TS U e b ; T s e T. 

b 

The requirement that the action S\ should be invariant under the group T implies several 
relations between the p a 's. In general, the e a 's span reducible representations of the group 
T . Depending on the number of these representations, some of the coupling constants p a 
may be equal. If there are three irreducible representations of dimension 1, all the p a are 
different. If there is one representation of dimension 2 and one of dimension 1, as it is the 
case in the triangular lattice where T is C 3v , two coupling constants are equal: pi = p 2 . 
In this case, since the action is quadratic in the fields, the invariance under the discrete 
group T is enlarged to a continuous invariance group 0(2) generated by: 




(ei,e 2 ,e 3 ) -> (ei,e 2 ,e 3 ) | -sm9 cos6> |, (8) 
and 

/ 1 \ 

(ei,e 2 ,e 3 ) -> (ei,e 2 ,e 3 ) 10 0. (9) 

\0 l) 

The action Si is thus 0(3) <8> 0(2) invariant in this case. 

Finally, if there is only one representation of dimension 3, as it is the case in lattices 
with tetragonal symmetry, p 1 = p 2 = p 3 and Si is invariant under G = 0(3) <E> 0(3). To 
summarize, depending on the symmetry group of the lattice, <Si can be symmetric under 
0(3) ® 0(p) with p = 1, 2 or 3. 

Since any rotation matrix R G SO (3) is given by an orthonormal set of three vectors, 
we can gather the e a 's into a rotation matrix R: 

R(x) = (ei(x)e 2 (x)e 3 (x)), (10) 

and the action <Si can be written into a different, but equivalent, form: 

S 2 = ^Jd D xTr{P{dR~ l ){dR)), (11) 

where P is the diagonal matrix: P = diag(pi,p 2 ,p 3 ) and R G SO(3). 

Using R, the symmetry operations on the e a can be written in a compact form. The 
action S 2 is invariant under left 0(3) transformations R — > UR , U G 0(3) and right 
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transformations: R — > i?V where V belongs to the 0(p) group which commutes with 
matrix P. We find again that depending on the value of p, S 2 is invariant under the 
group G = 0(3) <S>0(p);p = 1, 2, 3. The right 0(p) invariance reflects the original discrete 
symmetry of the microscopic Hamiltonian since it mixes the e a while the left 0(3) is 
the usual rotational symmetry. The discrete symmetry group of the Hamiltonian acts, 
in the continuum limit, as the 0(p) group. This is an artefact of the continuum limit 
and has no dynamical consequences since the number of Goldstone modes is given by 
the breaking of the 0(3) spin rotation group only. Indeed the number n of Goldstone 
modes resulting from the symmetry breaking G — > H, where H is the subgroup of G 
which leaves the ground state invariant, is equal to the number of broken generators of 
G: n = dim Lie(G) — dim Lie(H). In our case, the isotropy subgroup H consists of 
the transformations of G leaving the orthonormal basis e° invariant or equivalently using 
Eq.(||) and the above U and V transformations: 

H : R -> VR V = R , (12) 

where V 6 0(p) C 0(3) and is determined once V is chosen (see Eq.(^0] for an explicit 
expression of V in a particular example). This particular subgroup H is called the diagonal 
group 0{p)diag of a subgroup 0(p) in 0(3) times the 0(p) of G acting on the right. The 
symmetry breaking patterns described by action (pj) and equivalently by action (O) are 
therefore G/H = 0(3) (g> 0(p)/0(p)di ag with p — 1,2,3 depending on matrix P. These 
are all the possible symmetry breaking schemes that can undergo a frustrated Heisenberg 
spin model. In any case, there are three Goldstone modes. Before ending this section, let 
us emphasize that in the non linear sigma model, the dynamical properties depend only 
on the geometry of the coset space G/H and not on the field parametrization. Contrary 
to the Landau-Ginzburg model, the vanishing dimension of the order parameter in two 
dimensions allows infinitely many different parametrizations of the theory, differing in 
(possibly non linear) field redefinitions. The choice of one particular parametrization 
may be useful in discussing symmetry properties or renormalization group equations but 
does not change the physics. In particular we have seen that the actions S\ and S 2 
are equivalent. There is another equivalent form of S\ and S 2 that we shall use for the 
discussion of the 0(3) <g> 0(2)/ 0(2) case: 

S 3 = X - J d D x (<?! (Ve? + Vej) + g 2 ( ei Ve 2 - e 2 V ei ) 2 ) , (13) 

with: 

e a .e b = S ab , g 2 = ~-p 2 , g\ = Pi + Pi- (14) 

The latter expression of the action is obtained from (Q) by integrating over the constraint 
e3 = ei A e 2 and the relations among the coupling constants are derived in Appendix B. 

2.2 The particular case of the antiferromagnetic triangular lat- 
tice 

We shall consider the antiferromagnetic triangular lattice with iV-component spins, iV > 
3, as an example. Let us start by the case N = 3. The symmetry group of the system 
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is the product of the rotation group 0(3) acting on the spin components times the space 
group of the triangular lattice. The Hamiltonian density must then be G = 0(3) <E> C^ v 
invariant. The three spins Si, % = 1, 2, 3 of the elementary plaquette are co-planar in the 
ground state with the well-known 120 degrees structure. Then, we expect that only two 
vectors of the triad (ei,e2,e3) are necessary for the decomposition of the Sj. 

i) E = Si + S2 + S3 spans the trivial representation of C% v . This linear combination 
of the spins has a vanishing expectation value at T = 0. It cannot be an order parameter 
and corresponds to massive modes. 

ii) The two vectors: 



ei 



e 2 



/ __\/3+l 



oc 



Si 



2 



s 2 + s 3 \ 



V 



■y/3-1 
2 



Si 



2 



S2 + S; 



(15) 



/ 



span the two dimensional representation of C$ v . They have a non vanishing expecta- 
tion value in the low temperature phase and can thus be taken as an appropriate order 
parameter. The "local rigidity" constraint of section (2.1) means here: 



< E(V) >= 0. 



(16) 



This allows fluctuations of the spins between cells but not within the cells. This constraint 
is consistent with the symmetry since E is a scalar for 03^. Once the constraint is imposed, 
it is straightforward to show that ei(x) and 62(2) are orthonormal by use of Sf = 1. The 
action for the triangular lattice is then: 

Sy = \j d D x ( Pl ((Vei(a;)) 2 + (Ve 2 (*)) 2 )) . (17) 
Dombre and Read have obtained, from a direct microscopic derivation [?]: 



Pi 



V3J_ 



Note that the original C^ v invariance has been enlarged in Eq.flT7|) to a 0(2) group given by 
(HH): O = 0(3) ® 0(2). We will see that this action is not stable under renormalization. 
The most general renormalizable action which is compatible with the symmetry 0(3) <g> 
0(2) is given by equation dl3|) . Its general form is stable under renormalization so that 
we shall work only with it in the following. 

It is easy to generalize the above action to the case where the fields e a (x) have N > 3 
components. The symmetry group G is in this case 0(N) ® 0(2). The ground states are 
given by eq.@. Let us choose one of them, for example: 



( \ 
1 

V01/ 



(19) 



The unbroken symmetry group H of the low-temperature phase is then the set of matrices 
leaving this configuration invariant: 



S 



0(N-2) 








cos# 
sin 9 



— sin 9 

cos 9 



( \ 
1 

V o l J 



cos 9 sin 5 
— sin 9 cos # 



^0 0\ 
1 

V o l / 



(20) 



This group H consists in two subgroups: a 0(N — 2) and a diagonal 0(2). Therefore 
the action Si describes the symmetry breaking pattern G — > H = O(N) (g) 0(2) — > 
0(iV - 2) ® 0(2)^. 



3 Group theoretical construction of the non linear 
sigma model 

In the last section, we have derived three equivalent forms of the relevant action for 
Heisenberg frustrated magnets Si, S 2 , S3 (eq.(|3] ]Il|Jl3D ). Each of them has its own interests 
and shortcomings. Action Si is closely related to the microscopic Hamiltonian while action 
S 2 is suited to the discussion of the symmetry properties. Finally action S3 offers a possible 
large N, N > 3, generalization which we shall discuss in detail. The particular form of the 
action is irrelevant near D = 2 since the RG properties depend only on the geometry of the 
manifold G/H [?]. These intrinsic properties will be formulated in the language of group 
theory which provides an abstract but powerful framework [?]. We shall first deal with the 
standard 0{N) /0(N - 1) model. The construction of the 0(N) ® 0(2) fO(N - 2) ® 0(2) 
model is presented after this warm-up example. 

3.1 The 0{N)/0{N - 1) partition function 

The partition function of the 0(N)/0(N — 1) model is [?, ?]: 




The functional delta selects the configurations of S(x) with unit length. We can take 
advantage of this delta to integrate out one degree of freedom in S(x). Let us choose u, 
u J = 1, collinear to the magnetization and write S( 

S(x) = a(x)u + 7r(x) ; 7r(x)_Lu ; a 2 + tt 2 = 1. (22) 

After integrating out cr(x), Z can be rewritten as: 

Z = I^D* exp (--L / d D x ((dvr) 2 + (cVT^) 2 )) . (23) 

The low temperature T perturbative calculation of ( p3| ) starts from small fluctuations 
around the ground state: < S >= u. They correspond to the excitations of the 7r-field 
and are the usual spin waves. The 7r's consist of the N — 1 Goldstone modes coming from 
the breaking of 0(N) down to the rotation group 0(N — 1) that leaves the ground state 
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invariant, i.e. the 0(N — 1) around the u-direction. Once the symmetry breaking pattern 
O(N) — > 0(N — 1) is given the NLcx model is entirely determined up to the coupling 
constant which in this case is the temperature. 

We present now a matrix formulation of the 0(N)/0(N — 1) model. Let us choose a 
ground state S° = u. We can write the S(x) field as: 

S(x) = R(x)S°, (24) 

where R(x) is the O(N) matrix sending S° onto S(x). The partition function Z can be 
rewritten as: 

Z = J DR exp (-— J d D xTr{K{dR- 1 dR))^j , (25) 

where 

K aP = S°S°. (26) 

In this last equation the indices are those of vectors (a, (3) = 1, .., N and R G O(N). The 
relationship between S(x) and R(x) is not bi-univoque since for any rotation matrix h(x) 
leaving S° invariant: 

hS° = S°, (27) 

one has: 

R(x)h(x)S° = R(x)S° = S(x). (28) 

As a consequence, the action is locally (i.e. gauge) right invariant under the transforma- 
tion: 

R h {x) = R{x)h{x) , heH = 0(N - 1). (29) 

Some degrees of freedom in R e O(N) are thus unphysical. To obtain a bi-univoque rep- 
resentation in terms of matrices we have to choose one unique element in each equivalence 
class R h , that is to fix the gauge. The set of these equivalence classes is the set of O(N) 
rotations up to a 0(N — 1) rotation: it is 0(N)/0(N — 1). We can easily find one element 
per equivalence class in terms of the physical 7r-field. Let us write R(x) and h(x) as: 

4 I) h{x)= ( k l) h'eO(N-l). 

The matrix A is (N — 1) X (N — 1), V and V are aiV- 1-component vectors and B is a 
scalar. We use relation (p9|) to eliminate as many degrees of freedom in R(x) as possible. 
It is convenient to choose: 

h\x) = A^VX 1 !. (30) 
This leads to exactly one element per class given by: 



We identify V by applying L to S , eq.(p4|): 



^ =£(;). (32) 



(7 
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The element L can thus be written entirely in terms of the ir fields: 



k.')^ 1 "^*'' (33) 

The set of L- matrices is such that to any it corresponds a unique L(tv) G 0(N)/0(N — 1). 
The quantity (L~ 1 dL) belongs to the Lie algebra Lie(G) of G = O(N) and we have: 

{L~ x dL) = (L~ 1 dL)c-H + (L~ 1 dL)n (34) 

where (L~ l dL)n is in Lie(H). The partition function Z can be finally written as: 

Z = [ Dn e-^ s , (35) 



S = ~\J d D x Tr{[{L~ l dL) G - H ] 2 ). (36) 

We have used the fact that if is a projector: K{L~ x dL) H = 0. The partition func- 
tion in Eq.(P5D is globally G-invariant and locally (i.e. gauge) if-invariant. Once a gauge 
choice is made (as in (^,^)) no if -transformations are allowed in (^,^) and the G- 
transformations are in general not compatible with the gauge choice, i.e. they do not 
preserve the form of matrices L. This means that a G-transformation must be accompa- 
nied by a if-gauge-restoring-transformation: 

L(n') = gL(n)h(g,n). (37) 

Thus, G is non linearly realized on the 7r-fields. This is completely different from the 
Landau- Ginzburg model where G is linearly realized on the 7Tj fields. 

Equation (55) is the general expression for the partition function of a NLct model 
defined on a coset space G/H. This coset space can be viewed as a metric manifold so 
that it is convenient to formulate the theory in the language of differential geometry. 

Since L~ x dL belongs to Lie(G), eq. (|34]) can be rewritten as: 

L- 1 d tl L = e I fl T I + u;T a , (38) 

where the T a 's are the generators of Lie(H) while the T/'s are generators in Lie(G) — 
Lie(H). and uj® are respectively the vielbein and the connection in the tangent space 
of G/H. Under ( p7|) they transform as: 

= h-\x)(elTj)h(x), 
= h-\x){^T a )h{x) + h- 1 (x)d fl h(x). 

The T/'s span a representation of H since: 

[T a ,Tj] = fai J Tj. (40) 
As a consequence, the e^'s span a linear representation of H. Using (pB]), action S in 




eq.(]36f) can be written as: 
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S= l -ld D x eiefaj, (41) 



2 

where T]u is the tangent space metric given by: 

VIJ = —Tr(KTjTj), (42) 

with K the projector on G — H. In the 0(N)/0(N — 1) case it is given by eq. fl26|) . In this 
case, the N—l generators T/s of Lie(0(N))—Lie(0(N—l)) span the vector representation 
of 0(N — 1) so that there is only one coupling constant: rju = rjSij. However, in 
general, T]u is a diagonal matrix with several different couplings. The number of these 
couplings is the number of quadratic invariants under transformation (|40|) constructed 
with the e^'s. For a symmetric space there is only one such invariant. This is the 
case of 0(N)/0(N — 1) for example. For a non-symmetric homogeneous space such as 
O(N) ® 0(2) /0(N — 2) <g> 0(2) this number is larger than one. The formula = e I i d fl ir' 1 
leads to the more conventional form eq. (|23|) of the action of the NLcr model defined on a 
coset space viewed as a metric space equipped with the metric gij(n) = e\e J -r\ij [?]: 

Z = Dn exp (~ / d D x g^d^dir^ . (43) 

Eq.(|41"D and eq.fl43|) provide alternative descriptions of NLa models defined on a coset 
space G/H in terms of purely local geometrical quantities of the manifold G/H such as 
for example the metric, Riemann and Ricci tensors. It is equivalent to work either on the 



manifold itself (43) or in the tangent space (41). For practical calculations, it is extremely 



convenient to use the tangent space formulation we have discussed above. It can be shown 
that the geometrical quantities such as the Riemann tensor depend only in tangent space 
on the Lie algebras Lie(G) and Lie(H). More precisely they depend on the structure 
constants defined by the following commutation rules: 

[T a , Ti] = fa/Tj, 

[T a ,T b ] = f ab c T Cl (44) 
[T I ,Tj}=f IJ K T K + f IJ a T a , 

where T a e Lie(H) and Tj G Lie(G) — Lie(H). 



3.2 The 0(N) <g> 0(2) /0(N - 2) <8> 0(2) diag partition function 

In the case of the 0(N) ® 0(2) /0(N — 2) <g> 0(2)di a g model, the order parameter is the 
set of iV-component vectors (ei,e 2 ) and the action is given by action S3 eq.(|T^). Let us 
define the order parameter as the rectangular matrix: 

$ = ( ei ,e 2 ). (45) 
The 0(N — 2) ® 0(2) transformations can be written: 

*$' =* r (x) *$ *i?(x), (46) 
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where R G 0(N — 2) and r G 0(2). The ground state Eq.(|i~9l) is invariant under the 
transformations: 



*$° = h 1 (xf^°H(x), (47) 

with hi G 0(2) and 

= ( k V ) ■ »• e 0(W - 2). (48) 

Thus, the matrices r(x) and -R(x) are defined up to the following local transformations: 

— ► r(x) h\{x) 




R(x) H(x) (49) 
In the low temperature phase we can rewrite $ in terms of the 2N — 3 Goldstone modes: 

•<*> = ( W (,)&RF ) ' (50) 

where 7r is a (N — 2) x 2 matrix and u(x) G 0(2). The 7rf , % = 1, . . . , N, a = 1, 2 transform 
as two independent vectors under 0(iV — 2) and as a vector under 0(2). u(x)y/l2 — t 7i7i 
represents one extra degree of freedom which is scalar under both 0(N — 2) and 0(2). One 
can use the gauge freedom (f49|) to go from a general element R(x) <g>r(x) of 0(N) <g> 0(2) 
to the unique element in the same gauge orbit L ® 1 2 : 



L(tt(x),u(x)) = ( 1 , \ t ^ f . Z — ; — J . (51) 

The matrix L thus parametrizes the coset space 0(N) ® 0(2) /0(N — 2) ®0(2) diag . Note 
that this matrix is the same as for the coset space 0(N)/0(N — 2). 

In fact this is not accidental and we will use the following property to simplify our 
study: very generally the coset spaces G®X/ H®X<nag (where X is the maximal subgroup 
of G commuting with H) and G/H are topologically equivalent. We can thus work directly 
with the coset G/H keeping in mind that we search for an action which has G <S> X as 
symmetry (isometry) group. 

The vielbein of G/H defined as in eq. (|3"£D decompose into two irreducible representa- 
tions under the action of H £g> X. X itself spans the adjoint representation of X and is a 
scalar under H. G — H — X is irreducible because H®X is maximal in G, stated otherwise 
G/H ®X is a symmetric space. Thus, the two projected matrices (L~ 1 dL)\ G _ H _ x and 
(L~ 1 dL)\x transform independently under the right action of the H <S> X group so that 
there are two independent couplings % and r] 2 . We are thus led to the action: 

S = -I J d D x far^mh-H-x + V2tr(L~ 1 dL)* x ) . (52) 

Denoting by / the indices of Lie(G) — Lie(H), and among them by a the indices of 
Lie(X), this action may be rewritten as: 

S=\jd D x efefaj. (53) 
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where the tangent space metric r\u is given by: 



Vu = -vMT'T") - - r h )5 Ict 5jf S br{T a T p ). (54) 

We recall that in our case T l G Lie{0{N)) - Lie{0{N - 2)), T a G Lie(0(2)) and 
T a G Lie(0(N — 2)) and that the corresponding algebra is given in (f44"l). Action ([52]) 
is completely equivalent to the action S3 we have obtained from the continuum limit ( JT3| ) . 
We prove it in appendix B and derive the relations between the couplings gi,g2 entering 
in <S 3 and 771,772: 771 = gi/2 ; r) 2 = gi + 2g 2 . 



4 Renormalization of the NLa model in D = 2 + e 
4.1 General case 

The renormalizability in D = 2 + e of NLcr models defined on coset spaces G/H was 
studied by D.H. Friedan[?]. The (3 function gives the evolution of the metric gijij) with 
the scale: 

% (55) 



9/ hJir 



At two loop order it is given by the following expression: 

Pij{9) = ~ e 9ij + Rij + 7^TR ipqr R jpqr + 0(T 2 ). (56) 

where: R^ and Rj pqr are the Ricci and the Riemann tensors of the manifold Gj H equipped 
with the metric gij. 

In principle, it is enough to compute R^ and Rijki from the metric g^ to obtain these 
recursion relations. In practice, these calculations are tedious and some formal algebraic 
work has to be done first. The trick is to get rid in the calculation of any dependence 
on the coordinates ir l by going from the manifold itself to its tangent space, eq.(f|l]). The 
crucial advantage is that in tangent space, the Riemann and Ricci tensors are functions 
only of the structure constants of Lie(G) and that the tangent space metric r\u is 
constant, see eq.(f>§), and involves only the coupling constants. In the vielbein basis, 

Pu, (57) 



eq. (|55| , |56|) becomes: 

dl 

PiM = -equ + Ru + \t R ipqr R jpqr + 0(T 2 ). (58) 

The matrix 777 j is given in eq.fl54]) and the Riemann tensor in tangent space can be 
expressed as: 

RlJKL = f Ij a faKL + 2^ IjM (f MKL + fLMK ~ IkLm) 



cM 

L.J 



+ ^ {f I KM + fMIK — Ik Ml) (/j M L + flJ M — f 1 

T {fjKM + J'MJK — fKMj) [fl M L + /i/ M — f M Li) ■ (^9) 
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The indices a and {J, J . . .} refer to H and and G — H respectively. G — H indices are 
raised and lowered by means of rj IJ and rju and repeated indices are summed over. In 
NLcr models the j3 function and its derivatives allows to compute the fixed point and the 
critical exponent v. Note that since the (3 function is a tensor, it does not depend on a 
particular choice of coordinates. As a consequence, the mere existence of a fixed point 
as well as the value of the exponent v do not depend on the representation spanned by 
the order parameter. The other renormalization group function which is needed to give 
a complete description of the critical behavior is the Callan-Symanzik 7-function. This 
function is determined by the field renormalization Z: 

7 = — ^- ■ (60) 
From this function follows the anomalous dimension rj: 

v = livl^V) ~ e ( 61 ) 

where rjl, rfe are the fixed point values of the coupling constants. 

The factor Z is given at one loop order by the Laplace-Beltrami operator acting on the 
coordinate 7r\ It can be shown that this is nothing but g^T^ where is the Christoffel 
connection on the metric manifold G/H: 



Zn k = n k + -g ij Y%. (62) 



1 

— 1 

e' 

Once again, it is simpler to compute g^T^- by working in tangent space. We find for any 
coset G ® X/H <S) X where X is the subgroup of G that commutes with H: 



= -- (y, T aTa) n k + fE T ^l (63) 

where {T4} and {T a } are generators of G and X and where rji, r] 2 are defined in equation 
(|52j). J^aTaTa and ^ a T a T a are Casimir operators of G and X. In general, a choice 
of coordinates is not stable under renormalization. Equations fl5"2"D and ( p3| ) show that 
a good coordinate system which renormalizes multiplicatively consists in the 7r fields 
together with the massive a modes. They build up a linear representation of G <S> X 
such that the 7r's are an eigenvector of the Casimir operators with an eigenvalue that 
depends on the representation. Therefore, the 7 -function and thus the critical exponent 
n depends on the representation r ofG®X spanned by the order parameter. In our case, 
the Casimir operators have to be taken in the vector representation of both the O(N) 
and the 0(2) groups. Their values are therefore respectively N — 1 and 1. To summarize, 
in NLcr models the existence of a fixed point depends only on the symmetry breaking 
pattern G/H and not on the representation spanned by the order parameter. However, 
the universality class is completely determined once the representation r of G spanned by 
the observable is known. This scheme is completely different from what happens in the 
4 — e expansion where even the (3 function, and thus the mere existence of a fixed point, 
does depend on the representation of G spanned by the order parameter. In the following, 
we apply these results to the O(N) ® 0(2) /0(N — 2) <g> 0(2) models. For reasons that 
will soon become clear, we shall distinguish between the N = 3 and N > 3 cases. 
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4.2 Results for N > 3 

Using Eq.(|58|) we obtain the following two loop recursion relations valid for any A^ > 3: 



( 3 ^ i AT 



1% w 
2 Vl 8 



4 , N Tio 

-4 + 3 ( 1 -y)il 



m 



+(3JV-8)- 



(64) 



di]2 



N 



24 



.1 



Defining T 12 = 1/771,2, we find that, apart from the trivial zero temperature line of fixed 
points: T\ = T 2 = with Ti/T 2 arbitrary there is one non trivial fixed point C^l with 
coordinates: 

N - 1 / 1 3N 2 - lOiV + 4 



T* 



(iV - 2) 2 V 
1 (AT — l) 2 



[N - 2) 3 



-e 2 + 0{e A ) 



1 5N 2 - 16iV + 4 



(65) 



e 2 +0(e 3 ) 



z 2(iV-2) 3 \ 2 (N-2) 3 

This fixed point has one direction of instability so that our model undergoes an ordinary 
second order phase transition with critical exponent v: 



e + 



1 6N 3 - 21N 2 + 32N - 12 
(N — 2) 3 (2A r — 3) 



e 2 + 0(e 3 



(66) 



In order to complete our discussion, we have to specify the representation r of 0(N)®0(2) 
spanned by the observable of the physical system under study. We are interested in the 
AFT model with A^-component spins. In this case, the order parameter transforms under 
the vector representation of both O(N) and 0(2), see Eq.(^6|). At one loop, it follows 
from Eq.(^) and Eq. ([53]) that the anomalous dimension r] is: 



V 



4.3 Results for N = 3 



3N 2 - ION + 9 
2(A^-2) 3 ' 



+ 0(e 2 



(67) 



Although both the recursion relations and the values of the exponents given in the pre- 
ceding section are still valid in the N = 3 case the symmetry properties are less obvious. 
In this case, we can take advantage of the different equivalent parametrizations of the 
action we have derived in section 2. The convenient parametrization is given by Eq.flTTD: 



S 2 = I y d D x Tr (P(9i?- 1 )(9i?) 



(6* 



where P is the diagonal matrix: P = diag(pi,p2,P3) and R G 5*0(3). In the 0(3) <g> 
0(2)/0(2) case we have P\= Pi^ Pz- The relationship between the p^s and the tangent 
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space couplings r]u is given in Appendix B. Using Eq . (p4|) we can deduce the two loop 
recursion relations for the couplings pi. At the fixed point we find p\ = p* 2 = pi and thus 
P* oc 1. It follows from the discussion given in section 2 that the action <S 2 is 0(3) ®0(3) 
symmetric at the fixed point: the symmetry has been dynamically enlarged at the fixed 
point. Since 0(3)®0(3)/0(3) ~ 0(4)/ 0(3) the critical behavior of the 0(3)®0(2)/0(2) 
NLcr model is given by that of the 0(4)/0(3) NLcr model. It is a new result to find such 
a 0(4) symmetry for a Heisenberg system. We stress that it is not trivial to identify 
such a symmetry using a different parametrization such as the one given in action S3 (see 
Eq. (|T3"D ) . In this case, the 0(4) symmetry is non- linearly realized on the fields ei, e 2 . The 



critical exponents v and 77 are given by Eqs. (p6| , |67|) with iV = 3. Although the critical 
exponent v is identical to that of the N = 4 vector model, it is not so simple to get 77. The 
order parameter is a matrix R(x) = (ei(x),e 2 (x),e 3 (x))) (Eq.flTTD) and spans the tensor 
representation of 0(4). This point was previously missed in ref.[?]. As a consequence, 
the exponent rj of the Heisenberg AFT model is the anomalous dimension of a composite 
operator of the N = 4 vector model. To see this, we need the relationship between the 
0(3) matrix R and a 0(4) unit vector. It can be shown that to any unit 4-component 
vector: 

* = (*„,*i) ; *g + (69) 

i 

there exists a matrix R of 0(3) with components: 

R tj = 2(%* j - -<y + 2e ijk ^ ^k + 2(9 2 - ^)<% (70) 

Therefore, the expectation values of the vectors < e^(x) >,i = 1,3 are obtained from 
those of the bilinear forms < (^i&j — \o~ij) >■ 

We thus find no new universality class for Heisenberg canted models but instead the 
general phenomenon of increased symmetry at the fixed point. These models belong to 
the standard N = 4 Wilson-Fisher universality class. In dimension D = 3, the exponent 
v is very accurately known [?]: v = 0.74. However, the anomalous dimension of the 
composite operator (^^j — jSij) is only known at the two loop order in e = 4 — D. Let 
us finally emphasize that the phenomenon of increased symmetry at the fixed point is 
particular to the N = 3 case in the 0(N) <g> 0(2)/0(N - 2) ® 0(2) NLcr models. For any 
N > 3, the phase transition belongs indeed to a universality class different from 0(N) 
but as one reaches the physical N = 3 case one falls in the well known 0(4) one. This 
conclude our analysis of the NLcr models associated to canted magnets. 

The well-known e = 4 — D expansion starting from the upper critical dimension of 
the appropriate Landau-Ginzburg- Wilson (LGW) action has been applied to helimagnets 
more than ten years ago by Garel and Pfeutyf?] and Bailin et al.[?]. More recently, 
renewed interest on this subject has been drawn by Kawamuraf?]. We shall, in the next 
section, present the results obtained from this expansion. 



5 The linear theory and the e = 4 — D expansion 

The LGW action can be obtained in the same spirit as the NLcr. Once the symmetry 
breaking pattern is known, in our case 0(N) <g> 0(2) /0(N — 2) ® 0(2), all we have to 
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do is to find the most general action which is O(N) ® 0(2) symmetric and which ground 
state is 0(N — 2) Cg> 0(2) invariant. Among all possible actions, one has to select those 
which are renormalizable in D = 4, i.e. to keep only terms up to order 4 in the fields and 
to order 2 in their derivatives. The LGW action does not possess the invariance under 
reparametrization of the NLcr model. Moreover, only linear transformations of 0(N) <g) 
0(2) are allowed since non-linear transformations involve higher powers of the fields and 
their derivatives than allowed by renormalizability. For this reason the whole LGW or 
Linear theory depends explicitly on the representation of 0(N) ® 0(2) spanned by the 
physical order parameter. This fact have dramatic consequences on the renormalizability 
of LGW theories as compared to their corresponding NLcr models. In order to build the 
LGW action, we shall start from the NLcr model. In the particular case of canted models, 
we have to choose the parametrization which spans a linear representation of 0(N)<S>0(2). 
In this case the partition function Z is given by: 

Z = J De 1 De 2 5(ei.e 2 ) <5( ei 2 - 1) 5(e 2 2 - l)e~ 53 , (71) 

S 3 = ~ J d D x (<?! (Ve 2 + Ve 2 ) + g 2 ( ei Ve 2 - e 2 V ei ) 2 ) . (72) 

The ej are iV-component vectors. The LGW action is now obtained in a standard way 
by relaxing the constraints in Eq . (|71~D and use of a potential: 

V{ei, e 2 ) = \rn 2 {e\ + e 2 ) + Ul {e\ + e 2 ) 2 + u 2 {e 1 A e 2 ) 2 (73) 
The LGW action for canted magnets reads now: 

Slow = \j ' d D x Q (Ve 2 + Ve 2 ) + V{e u e 2 )) (74) 

We have rescaled the fields in order to obtain the standard normalization for the gradient 
term and have omitted the current term (eiVe 2 — e 2 Vei) 2 since it is not renormalizable. 
Note that it is precisely this term which allowed the NLcr action S% to be 0(3)®0(3)/0(3) 
symmetric at the fixed point when iV = 3. 

The two loop recursion relations for the couplings u± and u 2 were first obtained by 
Bailin et al.[?] and Garel and Pfeuty[?] . We shall here only summarize their results. Let 
us comment the RG flow: 

i) there is a critical value of N depending on the dimension: N C (D) = 21.8 — 23. 4e + 
0(e 2 ), under which there is no fixed point. In this case the transition is expected to be 
first order. Let us emphasize that for e = 1 the second term in the e-expansion of -/V c (e) is 
not a small perturbation of the first one since it is —23.4. This is the signal that a precise 
determination of N C (D) needs some control of the e-expansion which, as it stands, can 
not be used directly for e of order 1. 

ii) for iV > N C (D) there are still two different regions in the portion of the (ui,u 2 ) 
parameter space where the potential is stable: see fig. 1. One is the second order region. 
It lies above the line L joining the origin to an unstable fixed point (called C_) and is 
the basin of attraction of the stable fixed point, called Cl- The Heisenberg fixed point 
H is unstable towards Cl- The other region lies between the line L and the stability line 
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S of the potential: U2 = —2u\. It is a region of runaway behaviour and is expected to 
correspond to a first order region. We note that, as N tends to infinity, L tends to S, as 
expected. 



6 Interpolating between D = 2 + e and D = 4 — e. 

There is clearly a mismatch between RG results obtained in D = 4 — e dimensions from 
the LGW model and in £> = 2 + e dimensions from the O(N) g> 0{2)/0(N - 2) <g> 0(2) 
sigma model. Even though these NLcr models predict for any iV a continuous transition 
in the neighborhood of dimension 2, there are models with either N < N c or which do 
not belong to the basin of attraction of Cl for which the transition is expected to be of 
first order at least near D — 4. This is very different from the ferromagnetic case where 
both the 0(N)/0(N — 1) NLcr model and the LGW model predict the same critical 
behaviour. However, when iV > N C (D), there exists a domain in the coupling constants 
space M where a second order transition is predicted in both models. These domains 
are respectively the basins of attraction of C^l in D = 2 + e and of Cl in D = 4 — e. 
The natural question is whether or not these two fixed points are the same in a given 
dimension D between 2 and 4. The 1/N expansion allows to answer this question, at least 
for N large enough, since this expansion is non perturbative in the dimension [?]. The 
critical exponent v has been calculated to the lowest non trivial order in a 1/N expansion 
of the Landau- Ginzburg action (|74])[?, ?] : 



that i>i/ N (D) coincides with i/ 4 _ e (iV) and z/ 2+e (iV) to lowest order in 1/N. The same type 
of expansion can be done on the other exponents with the same results. 

We may thus conclude as in the ferromagnetic case that, when the fixed point exists 
near D = 2 and near D = 4, we can follow it smoothly from D = 4 — e down to D = 2 + e. 
Therefore, in the whole space E ={(M) = coupling constants, D,N) there should exist 
a domain Z where the transition is of second order and which is governed by a unique 
fixed point Cl{N, D) = Cnl{N, D). In the complementary of Z the transition is expected 
to be of first order. On the boundary V of these two domains, the transition should be 
tricritical in the simplest hypothesis. 

The situation can be summarized in the plane (N, D) of number of components of 
the model and dimension. The 4 — e findings have shown that there is a universal curve 
N C (D) separating a first-order region and a second-order region. If one believes that the 
2 + e results survive perturbation theory then the neighborhood of D = 2 belongs to the 
second-order region for all N > 3. As a consequence, the line N C (D) intersects the N = 3 
axis somewhere between D = 2 and D = 4. This defines thus a critical dimension D c that 
we do not expect to be a simple number. Making the hypothesis that the RG calculations 
have captured all the relevant fixed points then there are two possibilities: 




(75) 




4 - D or e 



D — 2, we find 



(76) 
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i) The critical value is between D = 3 and D = 4. This implies that the physical case 
N = 3, D = 3 undergoes an 0(4) transition as shown in (4.3). This case is favored by 
present numerical studies [?]. 

ii) The critical value is between D = 2 and D — 3. Then the physical case is governed by 
a fluctuation-induced first-order transition. It cannot be excluded that D c = 2 in which 
case the perturbative analysis of the nonlinear sigma model is always irrelevant. 

We note in addition that there is an additional possibility namely D c = 3. We do not 
see any reason why this would be realized since this looks like an artificial fine-tuning. In 
this case one may speculate that a tricritical mean- field like behaviour is seen for N = 3, 
D = 3. 

These alternatives are consistent with all the RG results and do not require additional 
fixed points not seen in perturbation theory. In this picture the stable fixed point seen 
in 4 — e for N > N c can be followed smoothly by the large-N limit till D = 2 and then 
identified with the conventional 0(4) fixed point via the 2 + e calculation of (4.2-3) for 
N = 3. 

7 Conclusion 

We have shown in this article that the non linear sigma model provides a new approach 
to the analysis of the critical behavior of frustrated systems. The double expansion in T 
and in e = D - 2 of the 0(N) <g> 0(2) /0(N - 2) <g> 0(2) NLct model has been performed 
and a fixed point has been found in D = 2 + e for any N, which turns out to have a 
remarkable 0(4) symmetry for three component spins. Since for iV < 21 the transition is 
expected to be first order near D = 4, we conjecture that for any N G [3, 21], the nature 
of the phase transition changes between 2 and 4 dimensions and is tricritical at the border 
of the second and first order region. 

In the simplest hypothesis the (N, D)-plane is divided in two region: a first-order region 
containing N = 3 and D = 4 and a second-order region containing the D = 2 line (for any 
N), the whole large-N line (for all D) and also in the neighborhood of D=4 the N > 21.8 
points. In between lies a universal line N C (D) whose 4 — e expression was already known. 
This universal line intersects the N=3 axis for some unknown critical dimension D c . If 
3 < D c < 4 then the physical point N=3, D=3 is second order in the 0(4) universality 
class and its exponent rj is that of a tensor representation. If 2 < D c < 3 the physical point 
is first-order. It may happen that D c = 3 (although we see no reason why) in which case 
one could see tricritical behaviour. This phase diagram is in agreement with all known 
RG results. To decide the fate of the physical point requires clearly additional techniques. 
Present Monte-Carlo results [?] favor 2 < D c < 3 although more work is needed. Direct 
RG calculations in D=3 may also help to confirm the phase diagram[?]. 

We note that in the generic case O(N) <g>0(2)/0(N- 2) <g>0(2), N > 3, the symmetry 
is not enlarged at the fixed point. In this respect iV = 3 is exceptional: for other values 
of N, the fixed point does not belong to the 0(N) Wilson-Fisher family. It would be 
interesting of course to investigate the fate of the XY N=2 case since known helimagnets 
have significant anisotropies that lead to a non-Heisenberg behaviour. 
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APPENDIX A 

We have seen in section 6 that for N = 3 the LGW and NLcr models do not predict the 
same critical behaviour. No fixed point is found in D = 4 — e dimensions from the LGW 
action while the NLcr model predicts a fixed point in D = 2 + e with a 0(3) ® 0(3) ~ 0(4) 
symmetry at this point. Actually, though these results are perturbative, it is easy to see 
that a 0(4)-symmetric fixed point can not be obtained with the LGW action ([74] ) since 
no value of (ui, u 2 ) makes this action 0(4)-symmetric. The reason is that the rectangular 
matrix {fa, fa) represents 6 degrees of freedom on which acts 0(3) on the right and 0(2) 
on the left and that this 0(2) can not be enlarged to 0(3) with only these 2 fields. This 
0(3) symmetry can be realized only with at least 3 fields: (fa, fa, fa). 

This would correspond to the 9-dimensional representation of 0(4). To check whether 
the discrepancy between the two models can be eliminated by allowing the LGW model 
to reach the 0(3) ® 0(3) symmetry, we have built and studied the most general action 
invariant under 0(3) <8> 0(2) and involving 3 fields: 



H = \{dfaf + \{dfaf + \{dfaf + \ Ul (0? + fa 2 f + \u 2 (0?0| - {fa.faf) 



(77) 



+\u'{<Plf + \u A <Pl (0? + fa 2 ) - \u, {{fa.faf + {fa- fa) 2 ) 

(0i, fa) is a doublet of 0(2) and 03 a singlet. H is 0(3) (g> 0(3) invariant when : 

u 2 = u 3 ; U\ + u 2 = «4 ; ui = u' (78) 

A one-loop RG calculation shows that no attractive fixed point exists in D = 4 — e in this 
model. This is the proof that the root of the problem is not only a question of symmetry 
but also a problem of dimension, field content and renormalizability. More precisely, if 
we set «2,«4 and u' to the values eq.([78|) which make hamiltonian eq. (|77|) 0(3) <8> 0(3)- 



symmetric, we find that the remaining symmetry in the broken phase is 0{3)diag- Then, 
the NLcr model associated with this symmetry breaking scheme is 0(4)/0(3). This NLa 
model is unique since it depends only on the Goldstone modes, and then only on the Lie 
algebras of 0(4) and 0(3) and not on the representations of these groups. On the other 
hand, there are as many associated LGW models as there are representations of 0(4) (or 
at least as there are actions built with representations of 0(4) that can be broken down to 
0(3)). Surprisingly, the LGW action built with the 4-component vector representation of 
0(4) admits a fixed point in 4 — e dimensions (the usual Heisenberg fixed point) and that 
built with the 9-dimensional tensor representation {fa, fa, fa) admits no such fixed point, 
as we have seen above. Since in these two LGW models, the symmetry breaking scheme 
is the same and then the Goldstone modes are the same, it means that the difference 
between these models lies in the massive modes. It is not clear up to now whether these 
modes can indeed be physically relevant for the critical behaviour. At least perturbatively 
and near 2 dimensions, the NLa model does not take care of these modes. In our case, 
we have to deal with the vector and tensor representations of 0(4) which both allow to 
replace the constraints of the NLct model by potentials and which lead to two different 
results. The relevance of the massive modes in the critical behaviour is then directly 
related to the way one chooses to go from the microscopic Hamiltonian to the different 
continuous actions, linear or non linear. 
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APPENDIX B 

We give in this appendix the relationship between different parametrizations of the 
sigma model. 

From tangent space to the constraints (for any N): 

We parametrize the matrix L in Eq.flSTD as: 

(79) 



(A 0i 2 ) 



where A is a rectangular N x (N — 2) matrix, 0i and 02 are iV-component vectors. Since 
L is in O(N), they must satisfy: 



l A.A = l N - 2 , t A0 1 =M^ 2 = O , 0f = 0^ = l , 0i.0 2 = O (80) 

A} A + 01. '01 + 02-*02 = Ijv (81) 

We are interested in the action of the NLtx model so that we have to compute L~ 1 dL: 



*A90i l Ad(j) 2 
L- l dL = | -*Ad<l>x fa.dfc 
tAdfc -01.902 



52) 



The projection of L l dL onto Lie(G') — Lie(H') leads to two sets of vielbein that are not 
mixed under the H-transformations: 



G'-H'-X 



/ M90i l Ad(f>2 

-*A<90 a 
V ',4<90 2 



(83) 



and 






>l-#02 





>l-<90 2 





J4) 



The total action is the sum of the traces of the squares of these two matrices weighted 
by 771 and r] 2 . These coefficients are by definition those coming from the tangent space 
metric rjjj. 



rjiTr 
T] 2 Tr 



G'-H'-X 
2 



-277! ((<90i) 2 + (<90 2 ) 2 ) + 4r / i(0i.«90 2 ) 2 



-277; 



15) 



'2W1- 



The coupling constants 771, r] 2 are now easily related to g\ , g 2 defined in 



V2 



9i/2 
9i + 2Q2 



56) 



From the P-matrix to the constraints (N = 3) 

We now compute directly Tr{P{R~ 1 dR) 2 ) to obtain explicitly the action (|TT|) . Let us 
start with R e 0(3), parametrized as follows: 



R 



(ei e 2 ei A e 2 ) 



57) 
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where e\ and e 2 are two 3-component vectors such that e 2 = e 2 = 1 and ei.e2 = 0. The 
diagonal part of {R~ x dR) 2 is: 



V 



-(de 2 y 



-(de l ) 2 -(de l ) 2 + 2(de l .e 2 f 



(88) 



so that: 



Tr (PiR-'dR) 2 ) = -( Pl +p 2 ) ((de,) 2 + (de 2 ) 2 ) + f(d ei .e 2 - de 2 .e x ) 2 (89) 
We obtain the relations between (pi,P2) and (gi,g 2 ): 

(pi =01 + 2a,, (go) 

1 = -2^2- 
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FIGURE CAPTIONS 



Figure 1: 

The renormalization group flow in the neighborhood of four dimensions for the Landau- 
Ginzburg model fl74[) in the case N > Nc{D). On the u\ axis one finds the conventional 
0(2N) Wilson-Fisher fixed point which is unstable towards the fixed point Cl- By fol- 
lowing a smooth path in the (N, D) plane we find that Cl is the 0(4) fixed point when 
N=3 in the neighborhood of D=2. 
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